Itinerant Spin Excitations in SrFe 2 As 2 Measured by Inelastic Neutron Scattering 
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We report inelastic neutron scattering measurements of the magnetic excitations in SrFe2As2, the 
parent of a family of iron-based superconductors. The data extend throughout the Brillouin zone 
and up to energies of ~ 260 meV. The spectrum calculated from a J1-J2 model does not accurately 
describe our data, and we show that some of the qualitative features the model fails to describe 
are readily explained by calculations from a 5-band itinerant mean-field model. In particular, the 
high-energy part of the spectra recorded above Tn do not differ significantly from those at low 
temperature, which is explained by the itinerant model and which has implications for theories of 
electronic nematic and orbital ordering. 

PACS numbers: 74.25.Ha, 74.70.Xa, 75.30.Ds, 78.70.Nx 



I. INTRODUCTION 

The new iron-based superconductors show an in- 
triguing interplay between structure, magnetism, and 
superconductivity-^— . It is likely that the paradigm to 
understand these materials will be rather different from 
that used to approach the cuprate superconductors, since 
the parent compounds for iron-based superconductors are 
bad metals, rather than Mott insulators. Theoretical es- 
timates suggest that the electron-phonon interaction is 
not the primary component of the pairing interaction in 
the iron-based superconductors 4 , and most attention is 
now on magnetic pairing mechanisms^. Much theoretical 
and experimental effort has therefore been devoted to un- 
derstanding the magnetism in this family of materials 6 . 
A key issue is whether the magnetism is better described 
within a weak-coupling (itinerant) 7 or strong-coupling 
(localized) 8 - picture. 

One of the most direct ways to identify the charac- 
ter of the magnetism is to study the magnetic excitation 
spectrum by inelastic neutron scattering (INS). Previous 
INS measurements of the antiferromagnetic (AFM) 122- 
arsenides AFe2As2 (where X is an alkali-earth metal) 
have been interpreted in terms of linear spin- wave theory, 
predicated on local-moment J]— J2 models^—. These 
models include nearest- neighbor (Ji) and in-plane diag- 
onal next-nearest neighbor (J2) Fe-Fe exchange inter- 
actions. A particularly interesting result from this ap- 
proach when applied to CaFe2As2 is the very large dif- 
ference between J\ a and Jib, the two in- plane nearest- 
neighbor exchange parameters^. Various mechanisms 
have been proposed to explain this anisotropy, including 
electronic nematic ordering^ 2 -, orbital ordering™ and the 
crystal structure itself 14 . An important piece of infor- 



mation, lacking up to now, is whether the anisotropy is 
modified on warming above the combined magnetic and 
structural transition temperature (Tn, s ), i-e. how the 
magnetic spectrum is modified on the change of sym- 
metry between low-temperature orthorhombic and high- 
temperature tetragonal phases. Recent resistivity mea- 
surements of de-twinned samples of Ba(Fei_ a; Co a ;)2As2 
show that there exists a large electronic anisotropy that 
persists above Tn,^. Previous INS measurements of 
CaFe2As2 above Tn,^ probed excitations only up to 
~ 60mcV. Given the possible role of magnetic fluc- 
tuations in the origin of the superconductivity of iron 
pnictides, further data on the magnetic spectrum in the 
tetragonal phase is of great interest. 

The questions we set out to answer here are, firstly, 
are the magnetic interactions in SrFe2As2 anisotropic in 
the ordered state, as in other 122-arsenides, secondly, 
do the spin excitation spectra change significantly on 
warming above Tn, S) and thirdly, how robust is the local- 
moment description of the INS spectra? The INS re- 
sults we describe here probe the magnetic excitations 
throughout the Brillouin zone (BZ) over the energy range 
5 < E < 260 meV, below and above Tn, s . We find that 
within a localized model the magnetic exchange param- 
eters are strongly anisotropic below and above 2n )S , and 
that an itinerant model gives the better qualitative de- 
scription of the data. The form of the data is not signif- 
icantly altered on warming above Tn,s- 



II. EXPERIMENTAL DETAILS 



Single crystals of SrFe2As2 were grown by the self- 
flux metho d 17 ' 18 . The crystals are highly homogenous, 
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as verified by microprobe analysis, a well-defined mag- 
netic/structural transition temperature of Tn s = 192 K, 
and by the observation of quantum oscillations from sam- 
ples from the same growth^. 

On cooling through Tn jS j the crystal symmetry changes 
from IA/mmm to Fmmm, and an AFM structure devel- 
ops with propagation vector Qafm = (0-5, 0.5, 1) referred 
to the IA/mmm unit cell [wave vectors are given in recip- 
rocal lattice units (r.l.u.)]. As the magnetic dynamics in 
SrFe2As2 are relatively two-dimensional we will hence- 
forth give only in-plane wave vector components. Fur- 
ther, from now on we index wave vectors with respect to 
the square unit cell formed by the Fe atoms in the para- 
magnetic phase (dFc = ^Fc = 2.8 A). In this convention 
Qafm = (0.5,0). The unit cells and their corresponding 
Brillouin zones are shown in Fig. [T] 




FIG. 1. (color online), (a) Projection of the real space crystal 
structure onto the afc-plane, with As atoms omitted for clarity. 
The solid (black) square is the unit cell of the Fe sublattice, 
the dashed (blue) square is the orthorhombic unit cell (space 
group Fmmm), and the dotted (red) diamond is the tetrag- 
onal IA/mmm unit cell, (b) Corresponding Brillouin zones. 
Referred to the Fe sublattice the standard square-lattice sym- 
metry points are Y = (0, 0), X = (0.5, 0), and M = (0.5, 0.5), 
and X is the AFM ordering wavevector. The arrows labeled a 
and /? indicate the direction of 1-d cuts and 2-d slices shown 
in Fig. Hand Fig. [3] 



The INS experiments were performed on the MERLIN 
time-of-flight (ToF) chopper spectrometer at the ISIS 
facility 1 -. Twenty one single crystals were co-aligned to 
give a mosaic sample of mass 5.4 g, with a uniform mo- 
saic of 4° (full- width at half- maximum). Spectra were 
recorded at temperatures of 6K, 212 K (Tn, s + 20 K), 
and 300 K with neutrons of incident energy E\ = 50, 100, 
180, 300 and 450 meV. The sample was aligned with the c 
axis parallel to the incident neutron beam, and the a axis 
horizontal. Data from equivalent positions in reciprocal 
space were averaged to improve statistics. The scattering 
from a standard vanadium sample was used to normal- 
ize the spectra, 5(Q, E), and place them on an absolute 
intensity scale, with units mbsr -1 meV -1 f.u. , where 
lmb= 10 m and f.u. stands for 'formula unit' of 
SrFe 2 As 2 . 
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FIG. 2. (color online). Energy-wavevector slices through 
the INS datasets at (a) T = 6K, (b) T = 212 K, and (c) 
T = 300 K. Panels (d) and (e) show simulations of the spec- 
trum over the same Q and energy range using a local-moment 
spin wave model with Ji a / J lb , J lb < 0, and Ji a = J lb re- 
spectively. Panels (f) and (g) show the calculated x"(Q,E) 
for the AFM phase, and the paramagnetic phase respectively; 
both taken from ref . [27l and convoluted with the instrumental 
resolution, with the energy rescaled by a factor 0.85. Data are 
from the run with Ei = 450 meV, and the intensity of both 
data and simulations have been multiplied by E to improve 
clarity. The Fe form factor has been included in the simula- 
tions. Dashed lines in panels (d) and (e) indicate the simu- 
lated dispersion relations for the high energy parameters (see 
caption to Fig. [3]) and for J\ a = Jib respectively. Dashed lines 
in panels (f) and (g) are the loci of maximum intensity for the 
calculated x"(Q,-E)- (The calculated low energy incommen- 
surate behavior in (g) is discussed in ref. [27] and likely arises 
from limitations of the model. See text for further discussion 
of this point.) 



III. RESULTS AND ANALYSIS 

The general form of the scattering at T = 6 K <C 
Tn, s is illustrated by the energy-momentum slice in 
Fig. Ufa). The magnetic spectrum is similar to that of 
CaFe2As 2 10ll:L , with intensity dispersing out of the Qafm 
positions. The dispersion is revealed in more detail in 
Fig. [3l which presents wave vector scans at four different 
energies. A single peak centered on Qafm = (0.5, 0) at 
low energies develops into a pair of peaks at ~ 75 meV 
which continue to separate and broaden with increas- 
ing energy. These peaks converge on Q = (0.5,0.5) at 
E - 230 meV. Figs.^b) and (c) show that at T = 212 K 
> Tn, s and T — 300 K ^ Tn, s the spectrum remains 
very similar to that of the ordered state, which means 
that strong AFM correlations persist well into the para- 
magnetic phase. 

Fig. [4] provides a comparison of the spectra recorded 
at 6K, 212 K and 300 K for wave vectors near to (0.5,0) 
and (0.5,0.5). The signals at 212 K and 300 K are some- 
what broader than that at 6 K, but otherwise there are no 
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FIG. 3. (color online). Wavevector cuts through the data 
taken at T = 6K with incident neutron energies (c) Ei — 
180 meV; and (a), (b) and (d) Ei = 450 meV. Solid (red) 
and dashed (blue) lines are calculated from fits to the local- 
moment Jia — Ju — J2 model with the low-energy and high- 
energy best fit parameters, respectively. For lower energies 
best fits were obtained with S e s = 0.30 ± 0.01, SJi a = 
30.8 ± lmeV, SJ lb = -5 ± 5meV, SJ 2 = 21.7 ± 0.4 meV, 
and SJ C = 2.3 ±0.1 meV. This is in contrast to higher ener- 
gies, where we found S cS = 0.69 ±0.02, SJi a = 38.7±2meV, 
SJ lb = -5 ± 5meV, and SJ 2 = 27.3 ± 0.7 meV. SJ C cannot 
be determined from cuts taken above the maximum of the 
dispersion along (0, 0, L) of ~ 53 meV. Dotted (black) lines 
indicate the estimated non-magnetic scattering at Q = (1, 0). 
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FIG. 4. (color online). Energy cuts at constant in-plane Q 
through the data taken with Ei — 450 meV. The left-hand 
panels show data at T — 6K, with (a) Q = (0.5,0.15) and 
(d) Q = (0.5, 0.5). The middle panels (b) and (e) show data at 
T — 212 K at the same two wavevectors. The right-hand pan- 
els (c) and (f) show data at T — 300 K, also at the same two 
wavevectors. Solid (red) and dashed (blue) lines are calcu- 
lated from the low-energy and high-energy fit parameters, re- 
spectively. The dotted (black) lines indicate the background, 
which is zero here because the non-magnetic background sig- 
nal from Q = (1,0) has been subtracted from these data for 
clarity. The calculations for different temperatures differ only 
in the damping fitted - the exchange parameters are the same 
for all six cuts. 



marked differences associated with the change of crystal 
structure, even well above Tj^s- 

We first compare the data with the linear spin-wave 
spectrum calculated from a local-moment Heisenberg 
Hamiltonian. This model is the standard one used 
to interpret the spin excitations in the parent 122 
compounds- - —, and is described in detail in the ap- 
pendix. The key exchange parameters in this model are 
Jia, Ju and J2, which define the exchange along the a 
and b directions and the ab diagonal respectively, as well 
as inter-planar exchange J c 

To fit data of the type presented we incorporated the 
resolution of the spectrometer in (Q, _E)-space, since the 
widths of peaks in the scattering are often very strongly 
coupled to the instrumental resolution, which is itself a 
function of Ei. We used the Toby fit software 20 , which 
uses Monte-Carlo methods to account for the effects of 
the (Q, E) resolution of ToF spectrometers. We fitted 
the measured scattering cross-section, in the form of a 
set of one-dimensional cuts, over the entire energy range 
for which magnetic excitations were extant, convoluting 
the cross section with the spectrometer resolution. 

Neutron inelastic scattering measurements of SrFe2 AS2 
performed on a triple-axis spectrometer with relatively 
low energy transfers 2 ^ reveal a single gap of 6.5 me V in 
the excitation spectrum. There is no evidence for two 
gaps, so we set the two modes described by equation 



Q in the appendix to be degenerate, i.e. C = in 
equations (j4j and the in-plane and out-of-plane single- 
ion anisotropy terms are equal (K a t, = K c = K). The 
single-ion anisotropy terms are thus determined by the 
size of the gap A, which is given by 



= K 2 ± [J la + 2 J 2 + J C )K (1) 



Because the c axis was parallel to the incident neu- 
tron beam, the (0, 0, L) component of the spectrum was 
coupled to time-of- flight, and hence to excitation en- 
ergy. However, by using several E\ one can determine 
the dispersion along (0,0, L), and hence determine the 
inter-layer exchange J c . We found, as expected from 
the calculated response functions^, that at low temper- 
ature the cross-section is a maximum for odd-integer 
values of L, and minimum for even-integer values, and 
that SJ C — 2.3 ± 0.1 meV. For energies above the BZ 
boundary energy along (0, 0,L), for which K = 0, 
{16(2SJ 2 + S Ji a )SJ c }^ 2 ~ 53meV, the INS spectrum is 
almost independent of L. Thus the fits performed where 
only data above ~ 100 meV were considered could not 
constrain J c , in which case we fixed J c to the value de- 
termined from the low-energy fits. 
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IV. DISCUSSION 

Superficially, the spin-wave model provides a reason- 
able overall description of the low-temperature data, as 
illustrated in Fig. [^d). One robust outcome from the 
analysis was a large difference between J\ a and Jib (see 
caption to Fig. [3]), similar to the anisotropy found in 
CaFe2As2^i. Another clear finding was that the damping 
term in the spectrum is energy-dependent, with a steady 
increase at low energies followed by a rapid increase at 
about 80meV. 

However, the spin-wave model fails in two important 
respects. Firstly, good fits could only be achieved by fit- 
ting the lower energy (E < 100 meV) and higher energy 
(E > 100 meV) parts of the low-temperature INS spec- 
tra separately - the low-energy parameter set is unable to 
account for the existence of an appreciable signal above 
— 200 meV, as seen in Fig. EJb) and Fig. H|d), while 
the high-energy parameters predict that the spin-wave 
branches below ~ 150 meV are unresolved, inconsistent 
with the data in Figs. [3][a) and (c). Secondly, the high- 
temperature spectra are inconsistent with the tetragonal 
symmetry, which constrains J\ a — Jn, for T > Tn, s - This 
is illustrated in Fig.[2je), which shows that the spin- wave 
spectrum when J\ a = Ju, is gapless at (0.5,0.5). The ori- 
gin of this softening is frustration. When J\ a = Ju,, the 
nearest-neighbor interactions are frustrated. The AFM 
structure can then be regarded as two decoupled, inter- 
penetrating AFMs, each on a square lattice of dimen- 
sions y/2ap c x v^iFc, the real space and Brillouin zone 
of which are shown as dotted (red) lines in Fig. Q] The 
magnetic unit cell and corresponding magnetic Brillouin 
zone for the uncoupled AFMs are indicated by the dashed 
(blue) squares in Figs. [Ha) and (b). The wavevectors 
X = (0.5,0), and M = (0.5,0.5) are both magnetic zone 
centers for the uncoupled AFMs, and therefore equiv- 
alent by symmetry. For the Fe sublattice, however, X 
and M are not equivalent by symmetry, and M is on 
the magnetic zone boundary for an itinerant AFM with 
Qafm = (0.5,0). So the fact that the spin- wave energy 
goes to zero at (0.5,0.5) for the tetragonal structure is 
purely a property of the local-moment treatment of the 
magnetic interactions. 

These results raise the question as to whether there 
exists a mechanism that maintains an anisotropic ex- 
change coupling ( Ji a y^z J lb ) in the paramagnetic phase, 
or whether the local-moment model must be abandoned. 
One possible mechanism is electronic nematic order, 
which has been widely discussed in connection with 
anisotropic properties of several other classes of strongly- 
correlated electron systems^— . The existence of an 
electronic nematic state has recently been proposed to 
explain anisotropy in the in-plane resistivity^ and elas- 
tic properties^ of iron pnictides. In some cases the elec- 
tronic nematic order is predicted to persist above Tn jS 
[ref. [Hj], although by only a few degrees at most. Thus, 
any evidence of broken symmetry above Tn, s is unlikely 
to be due to nematic order directly, but could perhaps 



be ascribed to nematic-like fluctuations, as in refi2£, in 
which it is suggested that nematic fluctuations might per- 
sist up to room temperature. Our higher temperature 
measurements were performed at both 20 K above Tn 
and room temperature, and as can be seen in Fig. [2jb) 
and (c), and Fig. SI the neutron scattering spectra are 
very similar. Although one might expect some change 
in the spectra over such a wide temperature range, we 
cannot completely rule out the existence of nematic fluc- 
tuations. 

Orbital ordering has also been suggested to explain 
the anisotropy^. In this case it has been predicted 
that a spin-wave mode would exist at high energies at 
Q = (0.5,0.5). The same calculations also indicate that 
on warming the energy of this mode softens slightly, 
and the peak height of the dynamical response function 
S al3 (Q, E) decreases rapidly. The broadening and lack of 
appreciable softening of the mode at Q = (0.5, 0.5) agree 
well with our data. However, as shown in Fig.@|c) to (e), 
the peak intensity is essentially unchanged, at variance 
with the calculation. Note that in ref. 13| the orbital 
order parameter is treated as an Ising-like parameter, 
so it is possible that if one considered instead a case in 
which there was partial orbital polarization the form of 
S a P(Q,,E) vs. temperature may be altered. 

If the local-moment model cannot be reconciled with 
our data, as seems likely, then what about itinerant- 
electron models^^ 7 , or hybrid models that combine both 
local moments and itinerant electrons^? In order to cap- 
ture the key features of the data it is likely that a rel- 
atively detailed model will be required. Indeed, a com- 
parison of the calculated \"(Cl,E) from a minimal band 
model (ref. l2q ) shows that it does not provide a good 
description of our data, except for < 50 meV. 

However, a mean-field model based on a more real- 
istic 5-band structure^ 7 - appears to give quantitatively 
good agreement with some of the features observed here. 
We illustrate this in Figs. HJa) and (f), which show the 
low-temperature INS data together with the calculated 
dynamical susceptibility x"(Q>E)> convoluted with our 
instrumental resolution. We also over-plot the locus of 
maximum intensity from the itinerant model calculation. 
This curve is rather more structured than the smooth 
dispersion curves shown in Figs. Hfd) and (e), with sev- 
eral abrupt changes of gradient, which would of course be 
parameterized in a local-moment treatment by exchange 
parameters that changed with energy. The energy scale 
of the calculated x"{Q>E) nas been changed by a fac- 
tor ~ 0.85 compared to the published calculation. This 
rescaling is likely due to the fact that in a mean-field ap- 
proximation the energy scale typically needs to be renor- 
malized down due to correlation effects not included in 
the model. The calculations in ref. [27] were performed 
with a Coulomb interaction U = 1.2 eV and Hund cou- 
pling J = 0.22 eV, chosen to yield the observed ordered 
moment of 0.8 fiB- There has, however, been some de- 
bate as to the strength of the electron correlations, as 
characterized by U and J, in iron pnictides22r— . 
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Fine tuning of the model parameters, i.e. U and J, 
would no doubt improve the description of the data, but 
the present set seems to work reasonably well. Thus, the 
spin fluctuation spectrum of SrFe2As2 that we measured 
indicates that this 122-arsenide exhibits rather weak elec- 
tron correlations. In addition we note that the calcula- 
tions from ref. 27 show that x"(Q,.E) does not soften 
at Q = (0.5,0.5) in the paramagnetic phase (Fig. Hfg)), 
which is also in agreement with our data (Fig. [2f b)). 
However, as noted in the caption of Fig. [3J the itinerant 
model of ref. does have some shortcomings, in par- 
ticular the form of the scattering at low energies in the 
paramagnetic phase. The calculation yields a signal at 
an incommensurate wavevector, which is not seen in the 
data. This incommensurability most likely arises from 
partial nesting in the paramagnetic band structure used 
in the mean-field calculation. Thus in order to obtain 
a better agreement between calculations and our results 
one would need to use a more detailed, experimentally 
determined, band structure, and also perhaps go beyond 
the mean-field approximation. 

The key advantages of the itinerant model are thus 
that: (i) it results in a more structured signal, as ob- 
served, which the local-moment model can only explain 
with discontinuous exchange parameters; (ii) it gives an 
explanation, in the form of particle-hole excitations, for 
the energy-dependent damping; and (iii) one does not 
need to invoke further phenomenology to explain the ab- 
sence of a soft mode at Q = (0.5, 0.5) in the paramagnetic 
phase. 

As mentioned above, there have been proposed re- 
cently hybrid models that combine both local moments 
and itinerant electrons^. In the AFM phase the calcu- 
lations yield a high-energy mode at Q = (0.5, 0.5), as we 
observe. They also yield a rather more structured disper- 
sion relation than a straightforward local-moment treat- 
ment, which also qualitatively agrees with our measure- 
ments. It is not clear, however, how the spectra would 
change on warming above Tn ]S , which is a crucial dis- 
criminant between local-moment and itinerant models. 



CONCLUSIONS 



In conclusion, our analysis shows that although super- 
ficially a local-moment model can be used explain the 
nature of the spin fluctuations in SrFe2As2 , close exami- 
nation shows that it fails in several respects. In particular 
the data cannot be fitted with a single parameter set, and 
a soft mode does not appear at Q = (0.5, 0.5) on warming 
above Tn, s - On the other hand, an itinerant model ap- 
pears to be able to explain both of these features. Thus 
it is not necessary to invoke further symmetry breaking, 
such as electronic nematic or orbital order, to explain the 
lack of soft mode at Q = (0.5, 0.5). 
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APPENDIX: LOCAL-MOMENT SPIN WAVE 
ANALYSIS 

The local-moment Hamiltonian used in the linear spin- 
wave analysis is: 



H = J2 JjkSj ■ S k + J2{K c (Sf) J+ K ab (S: : 

<jk) 



(2) 



The first summation is over nearest-neighbor and next- 
nearest-neighbor pairs with each pair counted only once. 
The Jjk are exchange parameters J\ a , J\ b and J2, along 
the a and b directions and the ab diagonal respectively, as 
well as inter-planar exchange J c . We also include in-plane 
and out-of-plane single-ion anisotropy constants K ab and 
K c . Diagonalization of equation © leads to two non- 
degenerate branches with dispersion 9 



(3) 



where 



A Q = 2S'{Ji fc [cos(Q ■ b) — I] + J la + 2J 2 + J c } 
+S(3K ai + K c ) 

C = S(K ab - K e ) 
D Q = 2S{ J la cos(Q ■ a) + 2 J 2 cos(Q • a) cos(Q • b) 

+JcC0S (Q_£) } . (4 ) 

Note that in cq. |4] a and b are the basis vectors of 
the Fe square lattice, whereas in Ref. |9| the correspond- 
ing equations are given with respect to the Fmmm space 
group. The in-plane cell parameters for the Fmmm space 
group are twice those of the Fe square lattice, as shown 
in Fig. Ula). Note also that in Ref. |9j a factor of two was 
missed from the term containing J c in the equation for 

The response functions per SrFe2As2 formula unit for 
magnon creation, ^"^(Q, ui), which relate to the neutron 
scattering cross section^, are given by 9 - 



A 



Q 



C-D, 



Q 



S"{Q,u) = S, 



C-D: 



{n(w) + l}*[w-wi(Q)], 



A 



clf- 



Q 



Q {n( W ) + l}^- W2 (Q)], 
(5) 
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where S c s is the effective spin and n(ui) is the boson 
occupation number. Only the transverse correlations 
(yy and zz for SrFe2As2) contribute to the linear spin 
wave cross section. In linear spin wave theory S e s = S, 
but we keep them distinct here because in the analy- 
sis the natural independent parameters are the energies 
SJ\ a , S Jib, ...SK a i,, SK C and the effective spin S e s of the 



fluctuating moment in the local moment approximation. 
To account for the finite lifetimes of the excitations we 
replace the delta-functions in equations ([5]) by damped 
harmonic oscillator functions^. Eqs. (A2-4) were used 
to calculate the neutron scattering cross-section using 
Ref.i, Eq. (4). 
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..e. we replace {n(E) + 1}5[E - £u, 2 (Q)] with — 



-E/k B T 



4rEi, 2 (Q)/7r 
[B 2 -E^ 2 (Q)] 2 +4[ri3]2' 



in S afi (Cl,E). 



